Electronic Structure of Samarium Monopnictides and Monochalcogenides 
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The electronic structures of SmX (X=N, P, As, Sb, Bi, O, S, Se, Te, Po) compounds are calculated 
using the self-interaction corrected local-spin density approximation. The Sm ion is described with 
either five or six localized /-electrons while the remaining electrons form bands, and the total 
energies of these scenarios are compared. With five localized /-electrons a narrow /-band is formed 
in the vicinity of the Fermi level leading to an effective intermediate valence. This scenario is the 
ground state of all the pnictides as well as SmO. With six localized /-electrons, the chalcogenides 
are semiconductors, which is the ground state of SmS, SmSe and SmTe. Under compression the Sm 
chalcogenides undergo first order transitions with destabilization of the / states into the intermediate 
valence state, the bonding properties of which are well reproduced by the present theory. 

PACS numbers: 71.20.Eh, 71.28.+d, 71. 15. Mb, 71.15.Nc 



I. INTRODUCTION 



The valency of rare earth compounds continues to be 
a vivid research ar e&m££££&l&& Of particular inter- 
est are systems where valency is influenced by control- 
lable external parameters like pressure, temperature or 
alloying. Several systems with varying or fluctuating 
valencies are knowniifi One of the most studied sys- 
tems is SmSjiiiaiSii 3 " which at low temperature and zero 
pressure crystallizes in the NaCl structure with a semi- 
conducting behavior, in accordance with the presence 
of divalent Sm ions and a full S p-band. At a mod- 
erate pressure of ~ 0.65 GPa SmS reverts to a metal- 
lic phase with a significant volume collapse of 13. 5%^ 
however retaining the NaCl structure. At very low tem- 
perature, the high pressure phase in fact reveals a small 
pseudogap, of the order of ~ 7 meVi Photoemission ex- 
periments show distinctly different spectra for the two 
phases, which usually are interpreted on the basis of diva- 
lent / 6 ions in the ground state and mixed valent / 5 — f e 
ions in the high pressure metallic phasei 1 ^* 1 ^ Recent ex- 
periments also find indications of small amounts of / 5 
character in the ground stated The C\i elastic constant 
decreases with pressure^ and becomes negative in the 
high pressure phase, indicative of an intermediate va- 
lence systemi The pressure induced semiconductor to 
metal transformation has also been traced by optical re- 
flectivity studies^ Mossbauer measurements^ resistiv- 
ity measurements^ and inelastic X-ray scattering. 8 Va- 
lence transitions can also be brought about by alloying 
with other trivalent ions, such as Y, La, Ce or Gd, into 
the SmS lattice* 1 *^ Similar valence instabilities are ob- 
served in SmSe and SmTe^iiiiiSiSli which also crystallize 
in the NaCl structure. For these compounds the volume 



changes continuously, but anomalously, with compression 
(at room temperature) i 2 ^ From the photoemission stud- 
ies it is concluded that SmSe and SmTe at amibient pres- 
sure, like SmS, are also of predominantly / 6 character^ 
while the monopnictides show clear signals of pure / 5 
ionsi 15 : 23 

The present paper addresses two issues of the SmX 
compounds. Firstly, density functional based total en- 
ergy calculations are used to investigate the ground state 
valency of Sm. To validate the approach, calculated val- 
ues of lattice constants and (for the Sm chalcogenides) 
valence transition pressures are compared to experiment. 
Secondly, electron spectral functions are calculated by 
combining the band structure Hamiltonian of the total 
energy calculation with an atomic description of the mul- 
tiplcts of the localized / electrons. 

The theoretical description of Sm compounds is a chal- 
lenge due to the Sm /-electrons. Conventional band 
structure calculations of SmSp* 2 ^ 5 -^ as implemented 
within density functional theory 2 ^ in the local density 
approximation (LDA) , describe the / electrons as narrow 
bands. A spin-orbit splitting of about 0.6 eV between oc- 
cupied 4/ 5 /2 and unoccupied 4/ 7 / 2 bands is found, and 
the valency transition can be traced by the crossing of 
the occupied / bands with the lowest d conduction band 
upon compressions 2 ^ However, the bonding of the / elec- 
trons when described as band states is greatly overesti- 
mated, as the calculated lattice constant of SmS from the 
LDA calculations is 7.6% smaller than the experimental 
valued To treat the / electrons as localized, Schumann 
et all 3 - applied self-interaction corrections (SIC) to the six 
4/5/2 states. This led to a semiconducting ground state 
of SmS, however with a too wide gap of charge transfer 
type, from S p to Sm d bands, in contrast to the Sm 
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/ to d character inferred from experiment A The reason 
is that the SIC scheme does not position the energy of 
the localized / states properly with respect to the band 
states. While the formalism does lead to some SIC / 
eigen-energies (which formally are just Lagrange multi- 
plier parameters of the theory with no physical interpre- 
tation), these are poor representations of the physical 
removal energies due to the localized nature of the exci- 
tations. To remedy this, subsequent work of Lehner et 
alii calculated the spectral density of SmS using a multi- 
band periodic Anderson model, in which the s, p, and 
d states were treated as band states within LDA, while 
the 4/ states were treated as localized states including 
atomic multiplets. The calculated spectral density was 
found to be in good agreement with the measured pho- 
toemission and inverse photoemission spectra. Recently, 
the electronic structure and optical spectra of SmS, SmSe 
and SmTe were calculated with the LSDA+U approach. 6 
This scheme includes an effective Coulomb parameter to 
separate bands of occupied and unoccupied correlated 
electrons. Some ambiguity exists as to how this is most 
appropriately done^ but the implementation of Ref. 6 
leads to good agreement with experimental optical spec- 
tra, although multiplet effects still have to be considered 
separately. 

In the present work the research efforts outlined above 
will be pursued further. The trends of the electronic 
structures of all of the samarium pnictides and chalco- 
genides will be investigated. We separate the issues of 
total energy and ground state determination from that 
of excited states properties. The total energy is calcu- 
lated using the self-interaction corrected local-spin den- 
sity (LSD) approximation. 29 When applied to rare earth 
systems this scheme may describe the electron states as 
either localized or delocalizedA For Sm compounds the 
relevant scenarios would have either 5 or 6 localized /- 
electrons on each Sm atom, corresponding to a trivalent 
or divalent Sm ion, respectively, and by comparison of 
the total energy the ground state configuration can be 
determined. The details of the SIC-LSD approach and 
its implementation will be presented in Section II. A of 
this paper. The SIC-LSD method is primarily suitable 
for calculating ground state configurations, crystal struc- 
ture and equilibrium lattice constants. Results for these 
parameters are presented in Section III, including the dis- 
cussion of pressure induced valency transitions. The pho- 
toemission spectra of SmX compounds are discussed by 
comparison to calculated spectral functions. The density 
of states from a total-energy calculation cannot be di- 
rectly compared with experimental photoemission spec- 
tra. This is particularly clear in systems with partly filled 
localized shells, where photoemission spectra reflect the 
multiplet structure of the final state. To include these 
effects, in the present work the LDA band Hamiltonian 
is combined with additional atomic multiplet informa- 
tion, which is calculated in an isolated reference atom 
in the Hubbard-I approximation^ The technical details 
of the present implementation 3 ! are discussed in Section 



II. B. In Section III the calculated spectral functions are 
presented and discussed in relation to experimental pho- 
toemission spectra. Finally, Section IV contains the con- 
clusions of the present work. 

II. THEORY 

Two important theoretical tools are employed in this 
work for the investigation of the electronic structure of 
SmX compounds. The ground state total energies are 
calculated with the self-interaction corrected local-spin- 
density method, 29,32 while the spectral functions are eval- 
uated with the 'LDA+-(-Hubbard F atomic corrections 
method. 31 

A. The SIC-LSD total-energy method 

The starting point is the total energy functional of the 
LSD approximation, which is renowned for its chemi- 
cal accuracy in describing conventional weakly correlated 
solids. 33 The essential approximation is the parametriza- 
tion of the contribution to the total energy, E xc , due to 
exchange and correlation effects among the electrons, by 
a simple term based on the homogeneous electron gas 
energetics. To facilitate an accurate description of the 
localized / electrons of rare earths, the self-interaction 
correction is introduced. This correction constitutes a 
negative energy contribution for an /-electron to local- 
ize, which then competes with the band formation en- 
ergy gained by the /-electron if allowed to delocalize and 
hybridize with the available conduction states. The SIC- 
LSD reduces the overbinding of the LSD approximation 
for narrow band states. Specifically, the SIC-LSD 2 ?. total 
energy functional is obtained from the LSD as: 

occ. 

E SIC = E LSD -Y^8 s a IC + E so , (I) 

a 

where a labels the occupied states and <5f IC is the self- 
interaction correction for state a. As usual, E LSD can 
be decomposed into a kinetic energy, T, a Hartree en- 
ergy, U, the interaction energy with the atomic ions, 
V ex t, and the exchange and correlation energy, E XC P- 
The self-interaction is defined as the sum of the Hartree 
interaction and the exchange-correlation energy for the 
charge density of state a: 

5 S J C = U[n a ] + E xc [n a }. (2) 

For itinerant states, <5f /c vanishes identically For lo- 
calized (atomic-like) states the self-interaction may be 
appreciable, and for the free atoms the SIC-LSD ap- 
proximation was demonstrated to be more accurate than 
LSD. 29 In Sm compounds, the self-interaction correction 
term is of the order S a ~ 80 mRy per /-electron. Fur- 
thermore, the spin-orbit term is added: 

E so = (^(f)t-s). (3) 
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We employ the atomic spheres approximation, whereby 
the crystal volume is divided into slightly overlapping 
atom-centered spheres of a total volume equal to the ac- 
tual volume. In (3), the angular momentum operator, 
I = f x p, is defined inside each atomic sphere, with r 
given as the position vector from the sphere center. Other 
relativistic effects are automatically included by solving 
the scalar-relativistic radial equation inside spheres. 

The E LSD functional includes the spin-polarization en- 
ergy, which is a very important contribution in the en- 
ergetics of SmX compounds. The Sm chalcogenides are 
non-magnetic^ while the pnictides are antiferromagnetic 
with very low Neel temperatures^*^ In the present work 
the spin-polarized energy functional is not used to de- 
scribe interatomic magnetism but rather to describe the 
intra-atomic exchange interactions. By Hund's first rule, 
the total spin of the / 5 or / 6 ion will be 5/2 or 3, and the 
ensuing energy gain is reasonably well described with the 
LSD exchange energy. In this spirit, the delocalization- 
localization transitions in elemental Ce^ and Am2£ were 
described in the LSD approximation as the points of on- 
set of spin polarization of / bands. In the SmX com- 
pounds the effects of a spin-polarized /-shell have negli- 
gible influence on the occupied non-/ band states, which 
are predominantly ligand p like. An important issue, dis- 
cussed recently for Eu monochalcogenides^ is whether 
ferro-magnetically ordered Eu / 7 ions can induce a spin 
splitting of the conduction band states, which can be ex- 
ploited in spin-filtering. 

Of Hund's second and third rules, the latter is governed 
by the spin-orbit interaction, Eq. (2), which also induces 
the formation of an orbital moment of the /-shell. The 
contribution to the total energy due to the second Hund's 
rule is, however, not fully accounted for by the spin-orbit 
interaction alone. The dominating contribution origi- 
nates from the energy gain by forming the eigenstates 
of the two-electron interaction 1/ry, the tetrad effect 
Little is known about how to implement these effects in 
self-consistent solid state calculations. The orbital po- 
larization scheme^ has had some success in describing 
the residual orbital moments of itinerant magnets like 
UFe 2 4i 

The advantage of the SIC-LSD energy functional in Eq. 
(1) is that different valency scenarios can be explored 
by assuming atomic configurations with different total 
numbers of localized states. In particular, these different 
scenarios constitute local minima of the same functional, 
E SIC in Eq. (I), and hence their total energies may 
be compared. The state with the lowest energy defines 
the ground state configuration. Note, that if no localized 
states are assumed, E SIC coincides with the conventional 
LSD functional, i.e., the Kohn-Sham minimum of the 
E LSD functional is also a local minimum of E SIC ' . The 
interesting question is, whether competing minima with a 
finite number of localized states exist. This is usually the 
case in /-electron systems^ and some 3d transition metal 
compounds,^ where the respective / and d orbitals are 
sufficiently confined in space to benefit appreciably from 



the SIC. The SIC-LSD energy functional in Eg. (1) is a 
true density functional, as discussed in Ref. \43L 

The SIC-LSD still considers the electronic structure 
of the solid to be built from individual electron states, 
but offers an alternative description of the single-electron 
states to the Bloch picture, namely in terms of periodic 
arrays of localized atom-centered states (i.e., the Heitlcr- 
London picture in terms of Wannier orbitals). Never- 
theless, there still exist states which will never benefit 
from the SIC. These states retain their itinerant char- 
acter of the Bloch form, and move in the effective LSD 
potential. The resulting many-electron wavefunction will 
consist of both localized and itinerant states. In con- 
trast to the LSD Kohn-Sham equations, the SIC elec- 
tron states, minimizing E SIC , experience different effec- 
tive potentials. This implies that to minimize E SIC , it is 
necessary to explicitly ensure the orthonormality of the 
one-electron wavefunctions by introducing a Lagrangian 
multipliers matrix. Furthermore, the total energy is not 
anymore invariant with respect to a unitary transforma- 
tion of the one-electron wavefunctions. Both of these as- 
pects make the energy minimization more demanding to 
accomplish than in the LSD case. The electron wavefunc- 
tions are expanded in linear- muffin-tin-orbital (LMTO) 
basis functions^ and the energy minimization problem 
becomes a non-linear optimization problem in the expan- 
sion coefficients. Further details of the present implemen- 
tation can be found in Ref. [32I 



B. Spectral Functions 

Photoemission spectroscopy is a powerful technique for 
investigations of rare earth compounds^ Several stud- 
ies of the electronic structure of Sm chalcogenides and 
pnictides by this technique have been reportedi 7 ^^*^*^ 
The photoemission technique probes the elementary elec- 
tronic excitations of the material, revealing energy- 
resolved and - through cross section variations - an- 
gular momentum decomposed information. The physi- 
cal excitation spectrum deviates substantially from the 
density of states (DOS) of the one-particle states ex- 
ploited for the minimization of the SIC-LSD total en- 
ergy. Strictly speaking, the DOS contains information 
on the one-particle band states that build up the SIC- 
LSD ground state of the SmX compounds. In view of 
Janak's theorem^ the LSD eigen-energies are good ap- 
proximations to total energy differences, provided the 
quasiparticle in question is of extended character. (An 
unknown shift exists, though, between addition and re- 
moval energies^ ) Therefore, for excitations which are 
extended, like the ideal Bloch waves of band theory, the 
LSD eigen-energies compare reasonably well with the 
physical excitation energies. However, often the photon 
impact creates localized excitations. In the SIC-LSD ap- 
proach one cannot expect the eigen-energies of the lo- 
calized /-electrons to have a straightforward physical in- 
terpretation. The true many-body character of the elec- 
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trons in the rare earth compounds is borne out in the 
multiplet structure of the photoemission spectra, which 
cannot be reproduced in an independent-particle picture. 
This shortcoming is also encountered with the LDA and 
LDA+U approaches. The multiplet structure arises from 
the two-electron interaction U ee — ^2[j 1/rij, where the 
prime implies a sum over pairs (i,j) of electrons. This 
interaction splits the ionic /" configurations into terms 
characterized by good quantum numbers of orbital, spin 
and total angular momentum, L, S and J. In the ground 
state of the solid usually only the lowest (Hund's rule 
coupled) term will be populated, maybe others as well 
if close in energy (within fegT, where fcs is Boltzmann's 
constant, or a typical / hybridization energy). But in 
the final state of the photoemission process, several of 
the / n_1 terms may be populated and contribute signif- 
icantly to the spectral features. 

Recently, Lichtenstein and KatsnelsorJ 1 devised a sim- 
ple procedure to augment LDA calculations to describe 
the multiplet features of partly filled shells. The main 
ingredient is an atomic selfenergy, T, atom (u)), extracted 
from the free ion Green's function: 

n atom ( \ _ {m\cn\n){n\c&\m) 



9n 



uj + E„ 



id 



(4) 



Here, m and n enumerate the ion multiplet states with 
energies E m and E n , c M and cj, are, respectively, the anni- 
hilation and creation operators of single /-electrons, and 
g mn is a weight factor specifying the relevance of m <-> n 
transitions. In a thermal environment, 



= -(e 
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(5) 



where (3 = l/fc^T and Z = J2 m e ~^ Em ^ s * ne partition 
function. The atomic levels \m) are found by diagonal- 
ization of U ee in the /" manifold: 

H atom =< f n. i^fn.j > +n(eQ _ ^ (6) 

where eo is the bare /-level and [If is the Fermi level. The 
matrix elements of H atom are expressed in terms of Slater 
integrals, F k , k — 0, 2, 4, 6 and Gaunt coefficients^ and 
the Hamiltonian matrix is easily diagonalized numeri- 
cally. The Slater integrals are calculated using the / 
partial waves of the self-consistent LDA calculation. 

From the atomic Green's function in Eq. (4) the 
atomic self-energy is extracted: 

and subsequently combined with the LDA Hamiltonian 
to give the solid state Green's function: 

1 



G solid i \ 
k ( LJ ) 



fjLDA 



(8) 



The spectral function is evaluated as the imaginary part 
of the Green's function of the solid: 

1. 
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(9) 



With care, the spectral function in Eq. (9) may be com- 
pared with experimental photoemission spectra. The lat- 
ter are also strongly influenced by matrix elements as well 
as secondary electron losses and shake-up effects, all of 
which are not considered in the present work. When com- 
bining the atomic and solid state electronic structure in 
this fashion, it is important to ensure that no interaction 
is counted twice. In particular, since the atomic multi- 
plets include the exchange interaction, the band Hamil- 
tonian in Eq. (8) should not be spin-polarized. Also, the 
LDA Hamiltonian includes the mean Hartrcc potential, 
which should therefore also be subtracted in Eq. (6), but 
this correction is indistinguishable from the Fermi-level 
adjustment necessary to embed the atom in the solid. 
This Fermi-level adjustment is the only fitting parame- 
ter of the scheme. If more than two /" configurations 
are compared, in practice also a rescaling of the first 
Slater integral, F°, is necessary to achieve results com- 
parable to experimental data. This is due to screening 
in the photoemission process, which in general renders 
bare Coulomb interaction too large (by a factor > 2). A 
completely analogous screening effect is considered in the 
LDA+U approach. 6 Lichtenstein and Katsnelson also in- 
voked screening effects of the higher Slater integrals for 
the interpretation of TmSe photoemission spectra^ but 
this was not considered in the present work. The present 
approach is similar in spirit to, but differs considerably 
in details from, the approach of Lehner et al. employed 
for SmS. 4 



III. RESULTS AND DISCUSSION 
A. Cohesive Properties 

The calculated total energies as functions of unit cell 
volume are shown in figures 1 and 2, for SmAs and SmS, 
respectively. In each figure two curves are drawn, cor- 
responding to the two cases of five or six localized /- 
electrons on each Sm ion. In SmAs the lowest energy 
is found for Sm(/ 5 ) ions, while for SmS the lowest en- 
ergy is found when assuming localized Sm(/ 6 ) ions. The 
energy differences between the / 5 and / 6 minima are 
AE(f 5 - / 6 ) = -74 mRy for SmAs and +15 mRy for 
SmS. The positions of the minima and the curvatures of 
the / 5 curve for SmAs and the / 6 curve for SmS yield the 
theoretical equilibrium lattice constants and bulk mod- 
ulii, which agree well with experimental values, cf. Table 
I. 

To investigate further the electronic structures of the 
ground states found for SmAs and SmS we show in Fig. 3 
the DOS for these two compounds. For both SmAs and 
SmS, one observes that the ligand s- and p-bands are 
completely filled. In SmAs, two narrow Sm /-resonances 
are situated just above the Fermi level. These are the / 
majority spin bands for the two / orbitals which have not 
been treated as localized. The minority / bands appear 
further above the Fermi level. In SmS, similarly, there is 
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one / majority orbital which is not localized and there- 
fore appears as a narrow spin-up / band, in this case ~ 3 
eV above the conduction band edge. As a consequence 
the density of states for SmS reveals a gap of 1.3 eV. 
This is similar to the gap found in Ref. (~ 2 eV) and 
in agreement with the fact that SmS by experimental ev- 
idence is a semiconductor. However, the gap in the DOS 
curve in Fig. 3(b) does not correspond to the physical gap 
of SmS, which is of / — ► d character and only ~ 0.15 eV 
in magnitude^ The DOS curves of Fig. 3 do not show 
the position of the localized states which, as discussed 
in Section II. B, are not well represented by their respec- 
tive SIC-LSD eigen-energies. Instead, the localized states 
can be calculated by total energy differences between the 
ground state and a state with one / electron removed. To 
achieve this a supercell consisting of four SmS formula 
units was considered, and the total energy calculated for 
either all Sm atoms in the / 6 configuration, or three Sm 
atoms in the f 6 configuration and one Sm atom in the 
/ 5 configuration. The missing electron in the latter case 
is artificially compensated by a uniform positive back- 
ground, in accord with the scheme often employed for 
charged impurities in semiconductors. 48 With this ap- 
proach the lowest / — > d transition is found to occur at 
0.58 eV below the conduction band edge in the DOS of 
Fig. 3(b). This is in considerably better agreement with 
the experimental band gap, albeit still too large, which 
could be an effect of the limited size of the supercell. 

In Fig. 3(c), the density of states of the high pressure 
(/ 5 ) phase of SmS is shown. Compared to the situation 
in Fig. 3(b) one more electron per formula unit is avail- 
able for band formation. At the same time, due to the 
reduced screening of the Sm nuclear charge, the unoc- 
cupied / bands move down in energy. The conduction 
states are being filled up to the position of the / reso- 
nance, which pins the Fermi level. The inset shows that 
the first of the two spin-orbit split majority / bands be- 
comes partly occupied. This means that each Sm ion is 
in a configuration of mixed / 5 and / 6 character, which 
we interpret as the SIC-LSD representation of an inter- 
mediate valence state, in accord with the conventional 
wisdom for the golden phase of SmSi The representa- 
tion of the intermediate valence state, as built from a 
hybridized band on top of an array of localized / 5 ions, 
may be a too simplistic description of the true quantum 
mechanical ground state of intermediate valence, yet due 
to the variational property of density functional theory it 
still leads to a good estimate of the binding energy. The 
lattice constant at the minimum of the SIC-LSD energy 
curve in Fig. 2, calculated for the high pressure phase, 
a = 5.69 A, is in good agreement with the lattice con- 
stant obtained for the SmS golden phase, a = 5.70 A>i 
or a = 5.65 Ai£ We note, however, that the calculated 
energy balance between the / 5 and / 6 phases in Fig. 2 
is at variance with the occurrence of the isostructural 
transition in SmS already at ~ 0.65 GPa^ The com- 
mon tangent construction applied to figure 2 would lead 
to a transition pressure of ~ 6.5 GPa, i. e. ten times 



higher. It is quite conceivable, though, that a more accu- 
rate description of the intermediate valence state would 
lead to higher binding energy for this phase, and hence 
to a lower transition pressure. With increasing compres- 
sion in the high pressure / 5 phase the / band gradually 
depopulates, at a rate of drif /d In V ~ 0.5, in accord with 
the experimental observation of increasing valency of Sm 
in the golden phase of SmS with pressured 

From plots similar to those in Figs. 1 and 2 we can de- 
duce the / 5 -/ 6 energy difference for all SmX compounds. 
Figure 4 illustrates the trends in this quantity. The cal- 
culations reveal a strong preference of the f 5 configu- 
ration in the early pnictides, with the energy difference 
of 135 mRy per formula unit in SmN. For the heavier 
pnictide ligands, the / 6 configuration becomes more and 
more advantageous, and for Bi it is only 6 mRy higher 
than the trivalent configuration. Moving to the chalco- 
genides, already in the Sm monoxide the / 6 configuration 
is found to be most favorable, by 6 mRy, and in SmS by 
15 mRy. Hence, the SIC-LSD total energy predicts a Sm 
valency transition between the pnictides and the chalco- 
genides. This is not in complete agreement with the ex- 
perimental picture, according to which the divalent and 
intermediate- valent states are almost degenerate in SmS, 
while SmO is trivalent and metallici 4 ^*^ Thus, it appears 
that the SIC-LSD total energy functional overestimates 
the tendency to form the divalent configuration of Sm, 
by approximately 15 mRy, in SmS. Assuming this error 
is similar for all SmX compounds, this would imply that 
the calculated energy balance curve in Fig. 4 lies too 
high by approximately 15 mRy. Therefore, the figure 
also includes an indication of how the energy difference 
behaves when a ~ 15 mRy correction is applied (dashed 
line). This switches the balance in favor of trivalcncy 
for SmO, which seems to be in better agreement with 
experiments, both with respect to the lattice constant 
and metallicity of SmOi-2*^ In the systematic study of 
the rare-earth metals and sulphides in Ref. a simi- 
lar uniform calibration (of 43 mRy) was applied to the 
trivalent-divalent energy difference of all the rare earths 
and rare earth sulphides, in order to account for the ex- 
perimentally observed valencies. The different size of the 
calibrating energy shift can to a large extent be traced to 
the neglect of spin-orbit coupling in Ref. 5, and the in- 
clusion of this effect in the present work. In conclusion, 
the SIC-LSD total energy functional predicts correctly 
the trends in trivalent-divalent energy difference through 
the samarium pnictides and chalcogenides, but fails on a 
quantitative scale of the order of 15 mRy. Such an er- 
ror is quite reasonable given that the functional does not 
contain any explicit contribution from the formation of 
atomic multiplets (the tetrad effected) , which would lead 
to larger energy-lowering for an f 5 ion than for an / 6 ion. 

The total energies as a function of volume in the Sm 
monopnictides and chalcogenides are used to study the 
basic ground state properties such as equilibirium lattice 
constant and bulk modulus. The calculated properties 
are in excellent agreement with the experimental values 
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which are given in Table I. All lattice constants agree 
within ~ 1% with experiment. Figure 5 compares the 
calculated lattice constants of the SmX compounds with 
experimental values and the equilibrium lattice constants 
of the competing valency configurations. 

In view of the intermediate valence interpretation of 
the high pressure f 5 scenario of SmS depicted in Fig. 
3(c), one may wonder if the occurrence of the sharp res- 
onance just above the Fermi level in SmAs in Fig. 3(a) 
does not also imply a slight admixture of / 6 into the 
ground state of this compound. Integrating the reso- 
nance band, one indeed finds ~ 0.11 occupancy of this 
band, i.e. SmAs is found to have a slightly mixed valent 
character. The admixture of / 6 character increases to 
0.17 in SmSb and 0.30 in SmBi. Neither the experimen- 
tal photoemission spectra of SmAs nor SmSb show dis- 
tinct features of f 6 character but on the other hand 
it is unclear whether a ~ 10 % admixture of / 6 character 
could be firmly excluded from the spectra. We are not 
aware of experiments on SmBi, but the present predic- 
tion is that a significant fraction of / 6 should be present 
in this compound. The proximity of the / 6 states is also 
probed by doping experiments starting from the pnictide 
and alloying with chalcogenidesi^SiS For SmAs, dop- 
ing with S or Se does not indicate valence instability for 
S concentrations up to 40 %££ or Se concentrations up to 
30 %j2A which suggests that in fact the extra electron of 
the chalcogen is not transferred into Sm / states, speak- 
ing against the presence of the unoccupied / band right 
at the Fermi level. On the other hand, when S is doped 
into SmSb the valency of Sm is seen to decrease already 
at low concentrations of It seems that more exper- 
imental as well as theoretical research, including alloy 
calculations, is needed to elucidate this issue further. 



B. Valence transitions 

As seen in figure 2 the divalent state, with six local- 
ized / electrons on each Sm ion, is the ground state of 
SmS, however, with compression the intermediate valent 
phase with five localized / electrons and some additional 
band / electrons becomes more favorable. From a com- 
mon tangent construction, a transition to the f 5 phase 
occurs, at a pressure of 0.1 GPa. This is accompanied by 
a volume reduction of 11.1%. Note, that this transition 
pressure corresponds to the total energy curves calibrated 
by the 15 mRy correction discussed in Section III. A. 

In the case of SmSe a similar transition is calculated 
to occur at a pressure of 3.3 GPa with a volume reduc- 
tion of 9.8%. Experimental evidence shows that SmSe 
undergoes a continuous transition in the pressure range 
of 2.6-4 GPa^i or 3-8 GPa£ The present theory can only 
describe a discontinuous transition. It has not been re- 
solved whether the continuous volume change is due to 
the experiments being conducted at room temperature 
or is an intrinsic property of the quantum state of SmSe. 
The calculated transition pressure and volume collapse 



are in reasonable agreement with the experimental data, 
cf. Table II, the experimental volume jump being esti- 
mated by extrapolation of the divalent pl^-relation over 
the anomalous region. Note that the two recent exper- 
iments quoted above disagree considerably with respect 
to the pressure range over which the transition occurs. 

Similarly to SmSe, experimentally SmTe also exhibits 
a continuous divalent to intermediate valent transition 
with pressure. The present theory finds a discontinuous 
transition occurring at 6.2 GPa with a volume collapse 
of 8.4%. Again, these values are in good quantitative 
agreement with the experimental data, cf. Table II. 



C. Photoemission 

The spectral functions of SmX compounds are calcu- 
lated as outlined in Section II. B. Figures 6 and 7 show 
the spectral functions for SmAs and SmS, respectively. 
The chemical potential of the reference ion is chosen such 
that the ground state is / 5 ( 6 if) for SmAs and f 6 ( 7 F) for 
SmS with an energy separation to the lowest n—1 excited 
levels (/ 4 ( 5 /) of 4.0 eV for SmAs and f 5 ( 6 H) of 0.8 eV 
for SmS, respectively), to coincide with the experimental 
values for these energies. The Slater integrals are almost 
equal for the two compounds, F k = 23.9, 10.6, 6.5 and 
4.6 eV, respectively, for k = 0, 2, 4, 6 (evaluated with the 
/-radial wave at an energy given by the center of gravity 
of the occupied /-partial density of states). However, for 
the direct Coulomb parameter, F° = U, a screened value 
of F° = 7.1 eV was adopted instead of the unscreened 
value quoted above. 

The SmAs spectral function in Fig. 6 shows the 
four distinct peaks corresponding to the / 4 ( 5 L), L = 
D,G,F,I, final states in the photoemission process. 
These states agree well with the three-peak structure ob- 
served by Ref. l23l at binding energies of approximately 
-10.0 eV, -8.2 eV and -6.0 eV (presuming that the 5 F 
emission is too weak to lead to a resolvable peak). In the 
positive frequency region one observes the "' F peak just 
above the Fermi level, in accord with the unoccupied ma- 
jority /-band in Fig. 3(a). The / 6 final states of S = 2, 
which correspond to the unoccupied spin down bands in 
Fig. 3(a), are situated further up in energy, however now 
with a considerable spread due to the many allowed mul- 
tiplets. The position of the corresponding levels in the 
reference atomic calculation are marked in the figure. 

The SmS spectral function is shown in Fig. 7. The 
spectrum is now characterized by the low binding energy 
three-peak structure, which is also observed by several 
experiments^ 1 ^^ at binding energies -0.8 eV, -1.5 eV 
and -4.0 eV, and which is attributed to the 6 H, G F and 
6 P final states^ The latter state coincides with the S 
p-band, as also found in the calculations. The results in 
Fig. 7 are similar to those obtained by Lehner et ali^ 
Recent experiments^ show traces of Sm / 5 emission in 
SmS photoemission experiments, possibly also present in 
older works m& It is unclear, whether this is due to small 
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impurity concentrations or implies a more complicated 
ground state already for the black phase of SmS. The 
present total energy calculations have found SmS to be 
a purely divalent system. It is well known that doping of 
SmS can lead to the intermediate valence phase, charac- 
terized by photocmission spectra of both the high and low 
binding energy type^£*2£ By carefully tuning the chemi- 
cal potential of the reference atom in the present theory 
we can indeed obtain a mixed spectrum, corresponding to 
a superposition of the spectral functions of Figs. 6 and 7, 
in good agreement with the spectra recorded for SmAs- 
SmS alloys^ The unoccupied states of SmS have been 
monitored with bremsstrahlung inverse spectroscopy^ 
The spectra reveal two broad structures, approximately 
4.5 and 9 eV above the Fermi level, which are in good 
agreement with the positions in Fig. 7 of the 8 S and e X 
features, respectively. 

IV. SUMMARY 

The cohesive properties of SmX compounds are well 
described by the local density approximation to den- 
sity functional theory provided the self-interaction cor- 
rection is applied to obtain an improved description of 
the atomic-like / electrons. The bonding properties 
arc quantitatively in agreement with experiment as evi- 
denced by accurate lattice constants for both the triva- 
lent pnictides and the divalent chalcogenides. Regard- 
ing the energy balance between the trivalent and diva- 
lent configurations of Sm in the studied solids (i.e. be- 
tween localized / 5 and f 6 configurations), the SIC-LSD 
approach seems to underestimate the bonding in the lo- 
calized f 5 configuration by 10-15 mRy, which can be con- 
sidered a minor error. However, for an accurate descrip- 
tion of the isostructural valence transitions induced by 
pressure it is a substantial inaccuracy. Correcting for 
this error we obtain good agreement with high pressure 
experimental results for SmS, SmSe and SmTe. The high 
pressure phase of the Sm chalcogenides is described in the 
SIC-LSD one-electron picture as an array of Sm / 5 ions 
with an additional partially occupied /-band, leading to 



a total / occupation between 5 and 6. For SmO, this 
is found to be the ground state. A small expansion of 
the SmO lattice, corresponding to an effective negative 
pressure, would lead to a transition to the divalent and 
semiconducting phase. This effect could be explored in 
SmO-SmS alloying experiments. 

The occurrence of multiplet effects in the photoemis- 
sion experiments of SmX compounds is direct evidence 
that the simple one-electron picture does not suffice to ac- 
count for all physical characteristics in these compounds. 
We demonstrated that the inclusion of local atomic corre- 
lation effects provides much improved spectral functions, 
as seen by the close correspondence between the calcu- 
lated main peaks and experimental photoemission spec- 
tra. A certain degree of fitting (of chemical potential 
and effective Coulomb interaction, U) goes into this pro- 
cedure, which hence cannot be considered as 'ab-initio' 
as the density functional based total energy calculations. 
Therefore, substantial further theoretical developments 
would be needed for a fully parameter free calculation of 
photoemission spectra. 

Concluding, this work has investigated the degree of 
intermediate valence in Sm pnictides and chalcogenides 
as manifested in three physical properties, namely the 
cohesive properties, the pressure characteristics and the 
photoelectron spectroscopies. 
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Compound 


Lattice constant ( 


a.u) 




Bulk modulus (GPa) 




Present 


Expt. a 


Present 


Expt. 


SmN 


9.46 


9.52 


131 


- 


SmP 


10.99 


10.88 


68.7 




SmAs 


11 IS 

11. lo 


111 fif 

11.10 


O f .0 


oQ.Z^o.O) , fo.O 


SmSb 


11.90 


11.84 


46.7 




SmBi 


12.14 


12.01 


41.3 




SmO 


9.41 


9.34 9 


106.6 




SmS 


11.25 


11.25 7 


53.4 


42(3) c , 50.3 d , 47.6(5.0) e 


SmSe 


11.70 


11.66 


53.4 


40(5) c 


SmTe 


12.43 


12.46 


37.6 


40(5) c 


SmPo 


12.64 


12.71 h 


33.4 





TABLE I: Calculated lattice constant and bulk modulii of SmX compounds in the NaCl structure. The Sm ions are in the 
calculated ground state configuration of / 5 for the pnictides and SmO, and of / 6 for the other chalcogenides. 
a : Reference except where reference to other work is given. 

52. 



Reference 
Reference 
Reference 
Reference 
Reference 
Reference 

Smo.532POr 

Reference 



at 110 K. 



Compound 


Pt (GPa) 

Theory Expt. 


Volume collapse (%) 
Theory Expt. 


SmS 

SmSe 

SmTe 


0.1 0.65°, 1.24 c 

3.3 ~ 4 a , 3.4\ 3 - 9 C ,2.6 - 4 d 

6.2 2 -8 a , 5.2\ 6 - 8 C , 4.6 - 7.5 d 


11.1 13. 5 a , 13. 8 C 
9.8 8 Q , ll c , 7 d 
8.4 9 c ,7 d 



TABLE II: Calculated isostructural transition pressures, P t (in GPa), and volume changes (in %), of Sm monochalcogenides. 
Experimentally, the transitions of SmSe and SmTe (at room temperature) are continuous, while SmS exhibits a discontinuous 
volume change. The calculated transition pressures include the 15 mRy calibration of the total energy calculated for the high 
pressure phase (see text). 
a : Referenced 

b : Insulator-metal transition of Reference I2CL 

c : Present author's estimates from figures of Reference and d : Reference IHiL The volume changes for SmSe and SmTe are 
obtained by extrapolation over the transition range. 
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FIG. 1: SIC-LSD total energy versus unit cell volume for SmAs. The two curves correspond to trivalent and divalent Sm 
ions, realized by localized f 5 (triangles) and / 6 (stars) configurations, respectively. The vertical bar on the V-axis marks the 
experimental equilibrium volume. 
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FIG. 2: SIC-LSD total energy versus unit cell volume for SmS. The two curves correspond to trivalent and divalent Sm ions, 
realized by localized f° (triangles) and / 6 (stars) configurations, respectively. The data presented do not include the 15 mRy 
calibration for the / 5 — / 6 energy difference (see text for discussion). The vertical bars on the 1^-axis marks the experimental 
volumes of the black and golden phases. 
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FIG. 3: SIC-LSD densties of states (DOS) for: a) SmAs with localized f Sm ions (a = 11.18 a.u.), b) SmS with localized f 
Sm ions (a = 11.25 a.u.), and c) high pressure phase of SmS with localized / 5 Sm ions (a = 10.75 a.u.). Energies are given 
relative to the Fermi level. The units of the DOS is states per eV and per formula unit. The full line is the total DOS, while 
the dashed line shows the DOS projected onto the Sm atom. The insets in a) and c) show a blowup of the region in the vicinity 
of the Fermi level. 
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FIG. 4: Trivalent-divalent energy difference (in mRy per formula unit) of Sm compounds (full line). A negative sign implies 
that the trivalent state is favored. The dashed line marks the corrected energy curve (see text). 
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FIG. 5: Comparison of experimental and theoretical lattice constants (in Angstrom) of SmX compounds. Experimental values 
(see Table I) are marked with solid circles, while lattice constants calculated assuming a divalent (trivalent) Sm configuration 
are marked with stars (diamonds). 
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FIG. 6: The calculated spectral function of SmAs at equilibrium volume, a — 5.91 A. The full curve shows the /-contribution 
and the dashed curve the non-/ contribution. The energy is given relative to the Fermi level. The main lines are characterized 
by their final state characteristics, either As s, p-bands, or / n±1 multiplet term. 
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FIG. 7: The calculated spectral function of SmS (black phase, a = 5.95 A). The full curve shows the /-contribution and the 
dashed curve the non-/ contribution. The energy is given relative to the Fermi level. The main lines are characterized by their 
final state characteristics, either S s, p-bands, or / n±1 multiplet term. 



